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^ . Abstract 



^ ■ In the previous papers we derived cosmological equations for models with Gurzadyan- 

^O , Xue dark energy, have performed their qualitative analysis and, particularly, have 
revealed a remarkable hidden invariance in the models with respect to the sepa- 

■»^ ■ ratrix Vt^ep in their phase portraits. Now, new analytic solutions for these models 

^^ , are obtained, showing additional symmetries at various curvatures. The likelihood 

^— V \ analysis with supernovae and radio galaxies data and their characteristics (age, 

■^ ' the deceleration parameter), again demonstrate the crucial role of the ^sep for all 

Oj. models, in spite of the diversity of both, the initial equations and their solutions. 
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1 Introduction 



The formula for dark energy, derived by Gurzadyan and Xue [1] represents a 
scaling between the speed of light c, gravitational constant G and scale factor 
a of the Universe 
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where h is Planck's constant, Lp is Planck's length. The formula is obtained via 
estimation of the energy of the relevant vacuum fluctuations and remarkably 
fits the SN data for the current value of the dark energy. 

Certain observational aspects of the GX-formula (1) are discussed in [2]. This 
formula can be interpreted as indication on possible fundamental constants 
variation with cosmic time. In particular, in order to keep interpretation of 
the vacuum energy as the cosmological term one has to suppose that the 
speed of light is proportional to the scale factor as studied in [3]. In other 
model gravitational constant changes with time (see similar model following 
from holographic approach [6]). 

Based on the scaling (1) we have proposed a set of cosmological models [3] 
and performed their general analysis [4] using phase portrait technique. Al- 
though cosmological equations describing each model are very different, some 
interesting similarities were found. 

In particular, qualitative analysis shows that there is a separatrix in the phase 
space of dynamical variables and it corresponds to particular solution of cos- 
mological equations with usual matter density parameter equal to 2/3 for 
almost all cases. In this way all solutions with Qm > ^sep contain Friedman- 
nian singularity in the past and have Friedmannian asymptotic in the future 
as well, while all solutions with Q^ < ^sep begin with some positive and large 
scale factor and zero density. 

In this Letter we find new analytic solutions of the models and compare them 
with observations on supernovae and radio galaxies. We then perform like- 
lihood analysis and provide best-fit values of the density parameter for all 
models and also calculate such important quantities as the age and decelera- 
tion parameters. 



2 GX models 



In general, with variation of the gravitational constant, speed of light and 
cosmological constant A with time we have [5] 
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where k = —1, 0, 1 is curvature parameter, fi is matter density, p is pressure, H 
is Hubble parameter, a is the scale factor of the Universe. Introducing instead 
of J^ derivatives with respect to scale factor ^, with p = we arrive to 



dfx ^a dG I dA I dc ,^ ^ »n ^ z.n 

da a Gda SttG da AnG cda 



Then, as in usual Friedmannian cosmology we introduce density parameter 



SttG 
where Hq is Hubble constant. 
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Model I. Neither the speed of light nor the gravitational constant vary with 
time, but A oc a^^. Cosmological equations are 
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If one neglects the source term in the continuity equation (that is possible for 
sufficiently large scale factor), this system describes the usual Friedmannian 
cosmology with negative curvature. 

There is analytic solution for density which reads 



where /x(ao) = /^o- Taking the limit of this expression with a — * one can find 
the value of Vim using (5) for separatrix [4] which turns out to be 



47r2 
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having values ^sep ~ 0.6345 for A; = — 1, Vts^p = 2/3 for A; = and ^sep ~ 
0.7022 for k = 1. Solution for the scale factor is 

±A2t = ^aAi{aAi + Bi) - y^aoAK^i + 5i) + 

,Di (Va^+VaoAi + BA (9) 



where 



Ai = c^{27r^ - k), Bi = —ao{2Gfioal - vrc^). (10) 

o 

This solution reduces to a simple form for separatrix (8) with the correspond- 
ing condition Bi = 



a = ao 



+ VAt. (Ill 



1/2 



Model II. Speed of light changes with time, c = [2:;^) a, the gravitational 
constant does not change and A=const. Cosmological equations reduce to the 
following system 



If one neglects the source in the last equation (that is again possible for large 
a), which consists of two terms, proportional to H^ and HA correspondingly, 
we come to the usual Friedmann-Lemaitre cosmology with cosmological con- 
stant A'. 

There is complete analytical solution for this model 



M«)=^o^ + 4^(l-^), (13) 



a(t) = -^ + ^±|^cosh(tv/5.)± 

(14) 
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where 



A2 = ^ao (47rG/io - A) , ^2 = A M - ^ j . (15) 

Separatrix for this model is given by 



^sep = o (16) 



independently on spatial curvature. 

Solution (14) has particularly simple form with Q^n = ^sep which is 



a = aoexp|M -— jAtl. (17) 



Model III. The gravitational constant changes, G = ^ ^ f^j , while the 
speed of light does not and A oc a~^. Cosmological equations are 

i,= + ^ = ?£!(£)Yl + J^). ^ + 3H, = 2hJi + ^).(IS) 

where hgx = g^^ =const. 

These are very different from the usual Friedmann equations. Even if we ne- 
glect the source terms in the continuity equation, the first cosmological equa- 
tion does not correspond to the first Friedmann one. 

Again there is solution for density 



fi = fXo- + 2fiGx(l--)- (19) 

a \ a J 

Notice similarity with model II. 

Solution for the scale factor looks exactly the same as (9) with 



27r2 
yl3 = Ai, ^3 = - ^ao{^io-2^iGx). (20) 

The separatrix is again given by (16) with i?3 = 0. We have already presented 
in [4] analytic solution for /u(t). 

Model IV. The speed of light varies as c = (^-^^^) a^^"^ and G'=const. 
Here A oc a^^, and cosmological equations are 



with P = ^ (GpGxf (l - l^)> PCX ^ ^^1^ =const. 



The source term in the continuity equation can be neglected for large a, but the 
first Friedmann equation contains the term a~^ which is again very different 
from the standard cosmology. 

There is complete solution for this model with 
and 

^ UQ zn oy j-i.4 -r UQU4 -\ 

where 



a = oo ± t\/ A4 + aoB4 + — r, (23) 



^4 = ^ao (^oao - /^) , B, = 2J^i2n^-k). (24) 



This solution reduces to (11) for i?4 = or, in other words, when Q^n -^ 1- 



3 Age 



From (2) in general case we have 



f da 

tn= — 
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a^ 



(25) 



where aj is either (for Q^n > ^sep) or the value of the scale factor for which 
the density vanishes (for Q^n < ^sep)- There are explicit solutions for models 
II and IV, so the age can be computed directly, in particular for model IV the 
age with Qm > ^sep is given by 



2 



^ ^^A, + aoB^-VA,j. (26) 

Coefficients A and B can be computed from (5) and A(ao) = 3iJQ(l — Qrn) 
with appropriate definitions of ficx and pcx- 

Results for all four models are represented in tab. 1 for selected values of Qm- 
We assumed Hq = 70h Km/(s Mpc). 
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Age of the Universe, deceleration parameter and maximal redshift (for Qm = 0.3) 
in models I- IV. il,sep is equal to 0.6345 when k = —1 for models I and IV and it 
is equal to 0.7022 when k = +1 for the same models. For all the rest cases it is 

^ ^sep — ■^/ '^* 

One can see from tab. 1 that the age is very different and is acceptable for 
most cases comparing to existing bound given by the age of globular clusters, 
in agreement with cosmic nucleo-chronology [7] 



^obs 



11.2 10^Yrs. 



(27) 



In particular, for models I and III the age is 14 GYrs quite independently on 
k for Qm = ^sep- The age is larger than (27) also for models II and IV for 
Qrn > ^sep- Notice also infinite age for separatrix in model II which follows 
from (17). No one of the models with Q^n = 0.3 pass the age constraint (27). 

Deceleration parameter characterizes the second derivative of the scale factor 
at present 



Qo 






■a(to). 
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Observational bounds on it can be found in [10] with 



- 1.35 < go < -0.15, 
within 3a CL. 



(29) 



In tab. 1 we provide values for this parameter in all four models under con- 
sideration. Notice that solutions for separatrix are "ballistic" with qo = and 
age 14 GYrs for models I and III and also the same for model IV with Q^n = 1- 
For models with de-Sitter asymptotic such as II and IV q^ is quite large and 
negative. Deceleration parameter comes to the region (29) only for model IV 
with Qsep < ^m < 1 with any spatial curvature parameter. 

In all four models the maximal redshift for Q^n = 0.3 is computed from 



1 + z 



Oo 



(30) 



and we find it Zn 



2.3 for models II-IV and a little bit higher for model I. 



Putting everything together, it is clear that model IV with Qm > ^sep passes 
all tests up to now. 



4 Supernovae and radio galaxies fit 



In what follows we use supernovae and radio galaxies data to put additional 
constraints on GX models. For this purpose we use the same data as in [8]. It 
consist of 71 supernovae from the Supernova Legacy Survey [9], 157 "Gold" 
supernovae of [10] and 20 radio galaxies of [11]. So in total we have 248 sources 
with redshifts between zero and 1.8. Dimensionless coordinate distance de- 
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Results of the fit. 



fined as usual (see, e.g. [12]) 



to 



Yia) 



Hotto f c{t)dt 



(31) 



and it can be expressed in a more convenient form as Y{z) with the help of 
(30). 

First of all, for standard ACDM model with Vim = 0.3, Vt\ = 0.7 we find 
X^ = 1.16 with 245 Dof. Then for each model cosmological equations are solved 
numerically and dimensionless coordinate distance Y{z) is calculated. The fit 
was obtained by the usual least square technique for all spatial curvatures with 
0.6 < Qm < 0.999 for all four models. Results of the fit for models I, III, IV are 
shown in tab. 2. The model II does not fit the data well and results for it are 
not included in the table. 

As one can see from tab. 2 for model IV best fit is obtained when Qm = 0.999 
and it is the same for different values of k. In models I, III the fit is better 
with smaller Q^- It is interesting to note that curves for models I and III with 
different k practically coincide with each other. This is due to the fact that in 
the cosmological equations the corresponding term is 1 — ^, and is close to 1 
with all curvatures. There is also another symmetry between model I, III. As 
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Figure 1. Best fit curves for models I, III and IV. The curves for all curvatures, 
A; = 0, +1, —1, do coincide. 



one can see from tab. 2 in case when k = solutions for Y{z) are the same 
for both models. Clearly, the best fit among models is given by model IV with 
high n^. 

In fig. 1 we show our results for models I, III, IV along with with 248 sources, 
including 228 supernovae and 20 radio galaxies. Since the character of solutions 
changes for Q^ < ^sep for models I, III we took flsep as a best fit. As already 
mentioned before, for models I, III our curves are the same so in fig. 1 they are 
shown as the same line. 



5 Conclusions 



Analytic solutions are obtained for all Gurzadyan-Xue models. Solutions for 
models II and IV are in explicit form. The most simple ones turn out the 
solutions for the separatrix Qgep for each model: power law for models I, III 
and IV and exponential one for model II. 

The presence of the separatrix, defining two families of solutions, appears 
equally crucial at comparison with observational data, for example, leading 
to the preference of higher matter density Qm > ^sep in all models. This is 
mainly, because the maximal redshift in models with Qm < Qsep is small, 
although the age can be acceptable, and on the contrary, low matter density 
models have shorter age. Note, however, that rather simple matter density 
models have been considered, mainly, to probe the value of the separatrix. 

Thus, the analysis confirmed the crucial role of the separatrix not only for the 
behavior of the solutions, but also at comparison with observations, which a 
priori is never trivial since the initial equations are quite different. 
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